z — )■ —z symmetry of spin-orbit coupling and weak localization in graphene 
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We show that the influence of spin-orbit (SO) coupling on the weak locahzation efi^ect for electrons 
in graphene depends on the lack or presence of z ^ —z symmetry in the system. While for z — >■ —z 
asymmetric SO coupling, disordered graphene should display a weak anti-localization behavior at 
lowest temperature, z — > —z symmetric coupling leads to an effective saturation of decoherence 
time which can be partially lifted by an in-plane magnetic field, thus, tending to restore the weak 
localization effect. 



The effect of spin-orbit (SO) coupling in graphene rep- 
resents an example of a stimulating theoretical study that 
remains difficult to detect experimentally. The form of 
intrinsic SO coupling in the graphene band structure, 
suggested by Kane and Mele has fuelled the the- 
ory of Z2 topological insulators, but its strength is too 
weak for detection by conventional spectroscopic meth- 
ods [ll-Q . Here, we show how the presence of SO coupling 
in disordered graphene could be manifested in quantum 
transport measurements. Specifically for graphene, the 
presence of SO coupling may not necessarily lead to anti- 
localization behavior known in semiconductors and met- 
als @, and this reflects the presence or lack of 2 — > — z 
symmetry in the source of SO coupling. 

In general, weak localization is very sensitive to sym- 
metry breaking in the electronic system or to scatter- 
ing involving electron spin since it is formed by elec- 
trons propagating along long diffusive trajectories [6nl6i|. 
The typical behavior of electrons in metals with strong 
SO coupling results in a pronounced weak antilocaliza- 
tion effect manifested by positive magnetoresistance, in 
contrast to simple metals and semiconductors where the 
weak localization magnetoresistance caused by the in- 
terference correction to conductivity is negative. For 
graphene with broken z — > — z symmetry, by a substrate 
or deposits, we find that, at the lowest temperatures, 
SO interaction leads to the conventional weak antilocal- 
ization. In contrast, for a z — — z symmetric system, 
SO coupling leads only to a saturation in the size of the 
weak localization correction rather than antilocalization, 
which can be taken for a saturation of t^{T) as T — >■ 0. 
Then, we analyze the influence of an in-plane magnetic 
field on the interference correction to conductivity for 
both z — — z symmetric and antisymmetric SO coupling 
scenarios and find that it lifts the aforementioned satu- 
ration of Tip . 

The breaking of z — > — z symmetry in pristine graphene 
is usually associated with the addition of a Bychkov- 
Rashba term Hbr [ll-fsllsl fT6l - [l9| to the intrinsic SO cou- 
pling fiKM dl, H, 0, [3 in the graphene Hamiltonian: 



H = v'Ep + hxM + hsR + U 

hKM — ^^zSz, hsR = 
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The last term in H accounts for Zeeman splitting due to 



describe three types of disorder on the honeycomb lat- 
tice: spin-independent perturbations, SO coupling with 
z — z symmetric perturbations, and z — >■ — z asymmet- 
ric perturbations, respectively. Here, we use a symmetry- 
based approach to determine how electronic spin may be 
combined with lattice and valley degrees of freedom. We 
focus on the region near the Fermi level which lies in the 
vicinity of two inequivalent corners of the Brillouin zone, 
known as valleys, with wave vectors K± = ±(47r/3a,0) 
where a is the lattice constant, and the momentum mea- 
sured from the center of a valley is p = fi,k — fiK±. The 
Hamiltonian ([l} operates in a space of eight-component 
Bloch functions $ =[0K+,A,ti <t'K+,B.t, 0K_,s,t: 0K_,yi,t' 
4>K+,A,i, (t>K+.B,i, 0K_,s,4,i <pK.-,A,i Consisting of two 
valleys K+/K_, two lattice sites A/B, and two spin 
com pon ents t/^l'■ We use three sets of Pauli matri- 
ces [Til [T2I to describe spin s — {sx,Sy,Sz), sublat- 
tice 'isospin' E = {T,x,T,y,'Ez) and valley 'pseudospin' 
A = {Ax, Ay, Az) [l^]. The matrices s,T,,A all change 
sign upon time inversion so that their products are time- 
inversion symmetric and T,aSj, SjAi may be used as a 
basis for a phenomenological description of static disor- 
der leading to SO scattering. 

The results of this symmetry-based approach are sum- 
marized in Table H] which shows how s, E, A may be com- 
bined to form irreducible representations of the planar 
group Cg'„ [21*1 which combines the point group Cgt, of 
strictly two-dimensional graphene with primitive trans- 
lations, as appropriate for the description of two val- 
leys K±. Matrices E and A are confined to the two- 
dimensional plane of graphene and their behavior under 
symmetry operations is impervious to the third spatial 
dimension. Thus, they are invariant under mirror reflec- 
tion in the graphene plane so that, in the absence of spin, 
they only appear in the representations that are z — ^ — z 
symmetric. The presence of spin, however, introduces a 
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TABLE I: Irreducible representations of the planar group Cg^ 
[2H . as provided by matrices Ea, Ai, and Sj. Representa- 
tions Ai, A2, Bi, B2, El, E2 are part of the point group 
of two-dimensional graphene Csd , representations E'^, E'2, G' 
incorporate primitive translations. 
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tivity. We assume that the Dirac-like Hamiltonian wS.p 
dominates the electronic behavior and that diagonal dis- 
order, luo(r) in Eq. ([5]), determines the elastic scattering 
rate, w Tq^ = ttju^/H, where 7 = p£/(27rfi^z)) is 
the density of states per spin, per valley [111, [3 ■ The 
current operator corresponding to the Dirac-like Hamil- 
tonian is momentum independent so that the current 
vertex entering the Drude conductivity is renormalized 
by vertex corrections. Then, the Drude conductivity is 
equal to cr = Ae^jD where the diffusion coefRcient is 
D — and the transport time is twice the scatter- 

ing time, Ttr = 2ro [1]. 

The weak localization correction Sa may be written in 
terms of disorder-averaged two-particle correlation func- 
tions known as Cooperon propagators Cl where index I 
refers to pseudospin (related to A describing valley de- 
grees of freedom), and s refers to spin (related to s). All 
the Cooperons that we consider are singlets with respect 
to sublattice isospin E because all isospin-triplet modes 
have relaxation gaps l/ro Then, Sa may be 

written in terms of a summation with respect to sixteen 
Cooperons consisting of combinations of spin and pseu- 
dospin singlet and triplets: 



pseudovector that lies in three-dimensional space: is 
even under z ^ —z reflection, but in-plane components 
Sx,Sy are odd. Thus, SO terms containing Sz appear 
in z ^ —z symmetric representations, terms containing 
Sx, Sy appear in z — ^ —z asymmetric representations. 

In Table m SzSz is an invariant of the point group of 
graphene representing intrinsic Kane-Mele SO coupling 
fiKM and T,xSx — ^ySy describes the Bychkov- 

Rashba term Hbr which assumes the exis- 

tence of a transverse field breaking z ^ —z symmetry. 
The entries in Table U take into account possible SO scat- 
tering mechanisms due to defects in graphene: V^ym in- 
cludes terms proportional to Sz, and 14sy includes Sx,Sy. 
The term U, Eq. ([2]), describes disorder decoupled from 
the spin degree of freedom: uo(r)I describing the influ- 
ence of remote charges, Uz^z{i^)'^z^z describing different 
on-site energies of the A/B sublattices, and Ux.z{^)^x-^z, 
Mj, 2(r)I]j,A^ accounting for fluctuations of A/B hopping. 

The other terms in U, Ua,x{''^)^a^x and Ua^y{r)'Sa^y for 
a = x,y,z, generate intervalley scattering. We assume 
that different types of disorder in the Hamiltonian are 
uncorrelated and x-y isotropic: 

{Ua,l{v)Ua'M {v')) = It^ ;(5aa"5n'(5 (r - r') , 
(Qaj (r)aa',i' (r')) = ^l^Saa'^jj^S {v - r') , 

(/3,,,(r)A,,,,(r')) = pyiv6,r5{v-v'). 

In the following study, we employ the standard dia- 
grammatic technique for disordered systems Q to cal- 
culate the weak localization correction 5a to the conduc- 
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Here, the factors cq = I, Cx = Cy = Cz = —1 take into 
account the fact that singlet and triplet Cooperons (of 
both spin and pseudospin) appear with opposite signs, 
and A is the vector potential of homogeneous external 
magnetic field, B = rotA (Bz = dxAy — dyAx). 

Inelastic dephasing is taken into account in Eq. ([5]) by 
and, in general, symmetry-breaking perturbations 
[such as those contained in the Hamiltonian Eq. (H])], 
contribute relaxation gaps T!, to the otherwise gap- 
less Cooperons C^, as quantified in terms of relaxation 
rates summarized in Table [Hi Then, the zero-field 
temperature-dependent correction, Sp{0), to the sheet 
resistance, where Sp (0) /p^ = —6a, may be written as 



Sp{0) 



e'p'D 



E 



CsCl 



{Dry 



-1/2 



2TTh ^ 

s.l—O^x^y^z 



CgCi In 



qdq 

-1 



Ci{q). (6) 



(7) 



and the magnetoresistance, Ap{Bz) = 5p{Bz) — Sp{Q), 



3 



TABLE II: Scattering rates, due to symmetry- breaking terms in the Hamiltonian Eq. ([T]), that produce relaxation gaps in 
the otherwise gapless Cooperons Cj where j refers to spin, / to vaUey. The relaxation rate of the intervalley Cooperons t^^ 
and the intervalley rate result from spin-independent disorder, intrinsic r^\f and Bychkov-Rashba r^^^ rates arise from 
coupling of spin and lattice, rates T~Jg, r^^oi t~\, t~\ account for coupling of valley and spin degrees of freedom. 
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B = 



where '0 is the digamma function. In Eq. ([S]) the in- 
fluence of an external field Bz is taken into account 
through discrete values = {n+1/2) / (Dtb) where tb = 
h/ (ADeBz). Used in conjunction with the Cooperon gaps 
listed in Table |lll Eqs. ([SJHl) provide a general de- 
scription of the weak-localization correction and corre- 
sponding low-field magnetoresistance in the presence of 
SO coupling in graphene, parameterized by six SO scat- 
tering rates r^jv/, Tgji, 'Tzv,cj '''iu.e' '^zv,oj '''iv,o- 

In order to analyze the influence of SO coupling in a 
realistic experimental situation [isl - fisj . we consider - in 
the rest of this paper - the spin-independent interval- 
ley scattering rate to exceed the decoherence rate 
and the rates due to SO coupling. This means that 
{r^,r^,r|} > and, thus, the valley-triplet Cooper- 
ons in Eqs. ([7][S]) may be neglected. Then, the six SO 
rates may be combined into just two relevant combina- 
tions: a rate t^I^ due to z —z symmetric SO coupling 
(terms Vsym and Hkm 
asymmetric coupling (T4sy and hsR) 



and a rate t„„1 due to z — ?• —z 
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Here, accounts for the Dyakonov-Perel [l^ spin re- 
laxation contribution, and the other terms for Elliott- 
Yafet [23] spin relaxation. 



The application of an in-plane magnetic field produces 
an interplay between SO coupling and Zeeman split- 
ting, as in semiconductor quantum dots [U, [25|. In- 
plane magnetic field introduces an additional term 
in the Hamiltonian 5H = (fta;/2)^s where w = 2^BB\\/h, 

ez — hid is the Zeeman energy and £ = {£x,iy, 0), \i\ = 1. 
This couples the spin-singlet Cg to the triplets C° and 
Cy. The spin part of the matrix equation for the valley 
singlet Cooperons (Cg, C°, C°, C°) = C° has the form 

/ n —iuiix ~iu}£y \ 

-iw4 n + r,-i 
-iuj£y n + T,-i 
n + 2T4/ 
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where the limit £2 ^ 00 of large Zeeman energy essen- 
tially means that 3> hr^^. 

In the absence of in-plane field, ez = 0, the low-field 
magnetoresistance Eqs. (|8l9p is given by 
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In the absence of SO coupling, Bso = ^asy — 0, Eq. (fT^ 
would describe negative magnetoresistance correspond- 
ing to weak localization [ill Il2j (lower dashed curves in 
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no spin-orbit 



FIG. 1: The low-field magnetoresistivity in the presence of 
z — >■ —z symmetric (left) or asymmetric (right) SO scattering, 
as compared to the absence of SO scattering (lower dashed 
curves). Solid curves show the influence of SO scattering, 



-1 

' sym 



= 25r^ and r^^y 



25t^ , respectively, with Zeeman 



energy = (top) to ez ^ t^^ (bottom). 



Fig. [IJ. In the presence of z — >■ —z symmetric SO cou- 
pling only, f^asy = 0, the contribution of the third term 
in Eq. (fT2|) is diminished, and the first and second terms 
cancel each other, leading to a suppression of magnetore- 
sistance for < ^Bso (upper solid curve on the left of 
Fig. [T]) vi^hich mimicks the effect of a saturated value of 



T,„^ + T^^,ly,. When z — >■ — z symmetry is 



ip ' ' sym ■ 

broken, Sasy 7^ and Bso 7^ 0, there is relaxation of all 
spin-triplets, and the second and third terms in Eq. (|12p 
are suppressed, leaving the first (singlet) term to deter- 
mine anti-localization behavior at T — >■ with positive 
magnetoresistance at low fields ^ Bg^y 
In the limit Cz ^ f>'T~J- , 



Ap=- 



2nh 



B, 



Fl 



B, 



+ Bs. 



(13) 



This result shows that for z — > — z symmetric SO cou- 
pling {Baay = 0), in-plane field partially restores weak 
localization at the lowest temperatures lifting the limi- 
tation of Tip discussed above. In contrast, for z — > — z 
asymmetric SO coupling, in-plane field changes weak 
anti-localization into a suppressed weak localization be- 
havior. The low-field magnetoresistance calculated using 
Eq. ([S]) for intermediate values of is plotted in Fig. [1] 
for z — — z symmetric (left) and asymmetric (right) SO 
scattering. 

To summarize, among the two extremes of SO cou- 
pling in graphene, z — > — z symmetric and z — > — z 
asymmetric, the manifestation of the latter in quantum 
transport resembles that observed in a 2D electron gas 
in GaAs/AlGaAs heterostructures, whereas the former 
is peculiar for graphene. Experimentally, the effect of 
z — >■ — z symmetric SO coupling can be taken for a de- 
coherence time "saturation" {t^{T ^> 0) — > Tgym) at low 
temperatures. Unlike inelastic decoherence, such satu- 
ration can be partially lifted by electron Zeeman split- 
ting induced by a strong in-plane magnetic field making 
the negative magnetoresistance Ap{Bz) sharper when 
Tj^^(T ^ 0) — ?> 0. It is necessary to mention that a 
similar behavior of weak localization magnetoresistance 
should be expected in magnetically contaminated con- 
ductors dl]. Spin-fiip scattering of electrons from local- 
ized spins leads to saturation of at the value of the 
spin-relaxation time whereas in-plane field freezes local 
moments thus suppressing spin-flip scattering of electrons 
and restoring the full size of the weak localization effect. 
However, the size of the in-plane field lifting the "satu- 
ration" of in these two cases is different: polarization 
of magnetic impurities requires /iB-B|| > kT whereas the 
suppression of the effect of z — — z symmetric SO cou- 
pling occurs when hbB\\ > hr^l^. 
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